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Abstract 

Bound and resonance states of helium atom have been investigated inside a quantum dot by 
using explicitly correlated Hylleraas type basis set within the framework of stabilization method. 
To be specific, precise energy eigenvalues of bound Isns (^5®) [n = 1 — 6] states and the resonance 
parameters i.e. positions and widths of states due to 2sns [n = 2 — 5] and 2pnp [n = 2 — 5] 
configuration of confined helium below N = 2 ionization threshold of He'^ have been estimated. 
The two-parameter (Depth and Width) finite oscillator potential is used to represent the confining 
potential representing the quantum dot. It has been explicitly demonstrated that electronic struc¬ 
ture properties become a sensitive function of the dot size. It is observed from the calculations of 
ionization potential that the stability of an impurity ion within quantum dot may be manipulated 
by varying the confinement parameters. A possibility of controlling the autoionization lifetime of 
doubly excited states of two-electron ions by tuning the width of the quantum cavity is also dis¬ 
cussed here. 

PACS numbers: 31.15.A-, 31.15.V-, 32.80.Zb, 37.10.Gh 
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1 Introduction 

The subject of atomic systems under spatial confinement is of immense interest among the researchers 
since the advent of quantum mechanics as the spectral characteristics of atomic systems placed under 
different confinements change appreciably compared with those of free atoms m- A number of 
phenomenological potentials have been used to model atoms within cavities [3], atoms under pressure 
[1], impurities in quantum dots or nano crystals [3], nanopores [QE], fullerenes [S] and foreign atoms 
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in liquid helium environment [9] etc. The study of quantum dots (QD) has got considerable attention 
in recent times due to its fundamental importance in theoretical researches as well as in fabricating 
new functional devices. The QD’s (or artificial atoms), in general, contains several electrons subjected 
to an external confining potential and they show similar structural properties as compared to pure 
atoms. The structural changes of the impurity atoms inside QD’s w.r.t. the parameters of confining 
potentials provide huge physical insight about the interaction of the atom with its surroundings. 
Although the bound states of confined hydrogen and helium atoms have been studied extensively by 
several researchers mw, very few attempts have so far been made towards the quasi bound or 
resonance states of one electron impurity atom in an isolated QD mm and also in case of confined 
two electron systems mm- Transformations of two-electron bound states to Feshbach and then to 
shape resonances depending upon a parameter of model rectangular well-type potential representing 
the QD have been studied by Bylicki et. al. [12]. Sajeev et. al. [T3] and Genkin et. al. m showed 
that the singly excited bound states of a two-electron atom become resonance states for appropriately 
chosen parameters of an external attractive spherical Gaussian type confining potential used to model 
the QD. 

In the present work, we have considered a spherically symmetric finite oscillator potential [^ IlSlfTB] of 
the type, 

K(r) = -Do(l + c^r)e-'’™^ (1) 


for modeling the QD confinement. Here Vq is the depth of the potential well and the cavity constant 
Cw is defined as, 
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where A is the width of the potential. By tuning the parameters Vq and A one can change the shape of 
the potential given by equation (1). Such type of two-parameter (Vq and A) potential provides much 
control and flexibility in modeling the size of a QD. When r ^ 0 i.e. near the center of QD , Vc{r) ~ 
and thus a harmonic nature is observed in the potential for a given cavity constant c^. But for large 
‘r’, it deviates from the harmonic behavior. In fact. The FO potential is quite similar in profile to 
that of Gaussian potential. At the same time, it facilitates the computation of matrix elements in 
a simple and efficient manner, especially when the Slater-type orbitals are used in constructing the 
wave function with appropriate boundary conditions for a confined system. This FO potential was 
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used by Winkler |15] to study the two-electron bound and resonant states of helium in QD where the 
electron correlation was not included initially in the optimized wave function. Even the inclusion of 
electron correlation could not remove the uncertainties in their calculations m- Later, Kimani et. 
al. |16j applied the restricted Hartree-Fock method to estimate the ground states of many-electron 
close-shell quantum dots modeled by the FO potential where the electron correlations were included 
approximately. Chakraborty and Ho |3] made a sophisticated approach to deal with this problem 
by expanding the wave function in single exponent Hylleraas type basis within the framework of 
stabilization method, but their work was restricted to only the lowest lying doubly excited resonance 
state 2s^ of helium. It is worthwhile to mention that an appropriate knowledge of resonance 

structure of few-electron QD with and without a central impurity atom will help to understand the 
electron transport phenomena occurring in real semiconductor QDs m- 

Under such circumstances, we have studied the resonance parameters of states originated from 
2sns and 2pnp [n = 2 — 5] conhgurations of QD confined helium below N = 2 ionization threshold 
of He'^ in the framework of stabilization method m by using explicitly correlated multi-exponent 
Hylleraas type basis set. This method was successfully employed by the present workers |18H22j for 
calculations of resonance parameters of different resonance states of the free and confined helium-like 
ions. In the present study, the resonance parameters of the states under consideration are estimated 
over a wide range of width (A) for a fixed depth (Vq) of the FO potential. The energy values of bound 
Isns states [n = 1 — 6] have also been been reported. Moreover, the positions of Is, 2s (^S) and 

2p i^P) states of He^ have been estimated for a comprehensive understanding about the structure 
of QD confined helium. The variation of ionization potential of QD confined He with respect to the 
width of the FO potential have been studied. It has also been shown that the potential given by 
equation (1) breaks the orbital angular momentum {1) degeneracy in Coulomb field for the energy 
levels of hydrogen-like atoms. Finally, we have shown that for a fixed cavity depth (Vb)) tbe widths of 
the resonance states show oscillatory behavior with respect to the width (A) of the quantum cavity. 
It has been noted that for higher excited states, such oscillations are more pronounced. The paper 
is arranged as follows: A brief discussion on the present methodology is given in Section H, followed 
by a discussion on the results in Section HI, and finally concluding in Section IV with a view towards 
further use of the present techniques in related studies of spatially confined atomic systems e.g. QD, 
presuure confinement, SCP confinement etc. 
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2 Method 


For any state of even parity arising from two electrons having same azimuthal quantum number, 
the variational equation [23] can be written as, 

2/1 i\/'ij-\2 




, dr IJ dr 2 J 

subject to the normalization condition 


dVr 


ri,r2,9i2 


= 0 


(3) 


fdVr 


ri,r2,9i2 


= 1 


(4) 


where the symbols used in equation (3) and equation (4) are the same as in reference [23]. The effective 
potential is given by, 


= j] 


2=1 


- 1 “ Vc (ri) 

Vi 


+ 


1 

ri2 


(5) 


The multi-exponent correlated wavefunction [20] considered in the present calculation is expressed as. 


firi,r2,ri2) = ^ ( 2 ) 


2=1 


^ ^ ^ Cimnr{ripri 2 + exchange 


l>0 m>0 n>Q 


+ EE 


* J 


gi{l)gj{2) Cimnr{r^ri 2 + exchange 


l>0 m>0 n>0 


( 6 ) 


where, 


ry,(i) = 


(7) 


where, ex’s are the non-linear parameters. Here, ri and r 2 are the radial co-ordinates of the electrons 
and ri 2 is the relative distance between them. In a multiexponent basis set, if there are p number of 
non-linear parameters, then the number of terms in the radially correlated basis is and, therefore. 


p(p+T ^ „ 

2 ^ y 


, where q is the 


the dimension of the full basis (N) including angular correlation will be 
number of terms involving ri 2 (Mj- For example, as we have used here nine non-linear parameters, 
the number of terms in the radially correlated basis is 45 and with 10 terms involving different powers 
of ri 2 , the dimension of the full basis {N) becomes 450. The values of the non-linear parameters are 
taken in a geometrical sequence: ai = (Tj_i 7 , 7 being the geometrical ratio |25j . The wavefunction can 
be squeezed or can be made more diffuse by changing the geometrical ratio ( 7 ) keeping cji constant 
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throughout. To have a preliminary guess about the initial and final values of nonlinear parameter cr, 
we optimize the energy eigenvalues of states below N = 1 ionization threshold of He~^ by using 
Nelder-Mead procedure [26]. The energy eigenroots are then obtained by solving the generalized 
eigenvalue equation, 

KC = EgC (8) 

where, ^ is the Hamiltonian matrix, 5 is the overlap matrix and E’s are the energy eigenroots. The 
wavefunction is normalized for each width (A) of the FO potential to account for the modified charge 
distribution inside the QD. Each energy eigenroot plotted against the geometrical ratio ( 7 ) produces 
the stabilization diagram. Subsequently, we can calculate the density of resonance states from the 
inverse of tangent at different points near the stabilization plateau in the neighborhood of avoided 
crossings for each energy eigenroot. The plot of calculated density of resonance states versus energy 
for each eigenroot is then fitted to a standard Lorentzian profile. The best fit, i.e., with the least chi 
square (y^) and the square of correlation near unity yields the desired position {Er) and width 
(T) of the resonance state. 

For each width (A) of the confining potential, the energy eigenvalues of ‘^S and states of the confined 
one-electron ion is obtained by using Ritz variational technique considering the wavefunction as, 

= (9) 

i 

Where, 7 ’s are the nonlinear parameters and C’s are the linear variational coefficients. For He'^{ns) 
states [n = 1 — 2], we have considered 14-parameter basis set whereas for He^{2p) state we have taken 
13 parameters in the basis. In both the cases, I is ranging from 0 to 4. All calculations are carried 
out in quadruple precision. Atomic units have been used throughout unless otherwise specified. 

3 Results and Discussions 

To construct the stabilization diagram corresponding to each width (A) of the FO potential, repeated 
diagonalization of the Hamiltonian matrix in the Hylleraas basis set of 450 parameters is performed 
in the present work for 400 different values of 7 ranging from 0.63 a.u. to 0.77 a.u. For a depth Vq = 
0.2 a.u. and a width A = 4.0 a.u. of the confining potential, a portion of the stabilization diagram 
for ^5® states of confined helium below N = 2 ionization threshold of He~^ is given in figure-1. It is 
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evident from figure 1 that there exist two classes of states: 

1. First few energy eigen-roots lying below He~^{ls) [—2.184879 a.u.] level are insensitive with the 
variation of 7 . This feature clearly suggests that these energy eigen-roots originating from Isns 
configurations of QD confined helium are bound i.e. stable against auto-ionization. 

2. Energy eigen-roots lying between He~^{ls) and He~^{2s) [—0.607849 a.u.] are sensitive with the 
variation in 7 and gives rise to flat plateau in the vicinity of avoided crossings of the energy 
eigenroots in the neighborhood of some particular energy values. This is a clear signature for 
the presence of resonance states of QD confined helium. 

The present calculated bound state energy eigenvalues {—E) of Isns (^5®) [n = 1 — 6] states of He as 
well as the He~^{ls) energies for different cavity widths (A) starting from a very low value of 0.001 
a.u. (corresponds to almost a free case) to a high value of 1000.0 a.u. are illustrated in Figure 2. It 
is to be noted that for very small cavity width A = 0.001, the Isns [n = 1 — 6] energy eigenvalues of 
helium and the He~^{ls) threshold energy are nearly identical to the corresponding energy eigenvalues 
of the free ions and they remain almost unaltered upto the cavity width A = 0.1 a.u. We can see 
from equation (1) that, for A ^ 0, c^, —)■ 00 and thus, I 4 —^ 0 which produce no effect of confinement. 
In between A = 0.1 a.u. and 10.0 a.u., the energy eigenvalues of helium decrease monotonically 
and ultimately saturates at {Eigns + 2 x Vq) a.u. In a similar fashion, the threshold energy He~^{ls) 
saturates at (Eis + Vq) a-u. This feature is physically consistent as we can note from equation (1) 
that for A — )■ 00 , the cavity constant c^, ^ 0, so that Vc{r) —)■ — Vq- Thus the one and two electron 
energy levels will undergo a downward shift by Vq and 2Vo respectively for A —)■ 00 . The variation 
of the ionization potential (in eV) i.e. the energy required to ionize one electron from the ground 
state (Is^) of helium atom is plotted against the width (A) of the cavity in Figure 3. In accordance 
with the variation of energy eigenvalues of helium and its one-electron subsystem, it is evident from 
figure 3 that, the IP is identical with the vacuum IP for low values of A while for high values of A, it 
increases by an amount Vq ~ 5.44 eV. It is thus evident from figures 2 and 3 that rate of variation of 
energy values of the ions are significant when the size of the confining cavity is of the order of atomic 
dimensions. It is also remarkable that the stability of an impurity atom can be controlled by suitably 
tuning the size of a QD i.e. the depth and width of the representing cavity. 

An enlarged view of the stabilization diagram (given in figure 1) for states of He within the energy 
range —0.8 a.u. to —0.64 a.u. is given in figure 4. The states of He below N = 2 ionization 
threshold of He^{2s) can arise due to 2sns and 2pn'p {n,n' > 2) configurations. From a closer look 
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at figure 4, we can see that for a short range of 7 , each eigenroot between N = 1 and N = 2 ionization 
thresholds of He'^ becomes almost flat in the vicinity of avoided crossings in the neighborhood of 
different energies. In order to calculate the exact resonance parameters, the density of states (DOS) 
p{E) is calculated by evaluating the inverse of the slope at a number of points near these flat plateaus 
of each energy eigenroot using the formula |18j given by: 


Pn{E) 


li+i - li-i 


En{li+l) 


En(pii—l) 


En{li)=Ei 


( 10 ) 


The estimated DOS Pn{E) is then fitted to the following Lorentzian form [18 
A ^ 


Pn(E) = yo + 


E(E-Erf+(^y 


( 11 ) 


where, yo is the baseline background, A is the total area under the curve from the baseline, Er gives 
the position of the center of the peak of the curve and represents the full width of the peak of the 
curve at half maxima. Among different fitting curves for each eigenroot corresponding to a particular 
resonance state, the best fitting curve i.e. with least and the square of correlation (R^) closer to 
unity [T 8 ] leads to the desired resonance energy (Er) and width (T). The evaluation of DOS following 
the fitting procedure has been repeated for each width of the confining potential (A). For example, 
the calculated DOS and the corresponding fitted Lorentzian for the 2s^ (^5®) resonance state of He 
below He~^(ls) threshold for cavity width A = 4.0 a.u. (given in figure 5) yields resonance position 
Er at —0.98163 a.u. and width T = 6.9961 x 10“^ a.u. 

The estimated resonance energies of doubly excited 2sns [n = 2 — 5] (^5®) and 2pnp [n = 2 — 5] (^5®) 
states of helium and corresponding 2s and 2p threshold energies for the cavity depth Vq = 0.2 a.u. 
and cavity width (A) ranging from 0.001 a.u. to 1000 a.u. are given in figure- 6 , while the variations 
of resonance energies (Er) of 2pnp [n = 2 — 5] yS^) states and corresponding 2s and 2p threshold 
energies versus A are given in figure-7. We have noted the following points. 


1. It is clear from figures -6 and 7 that for A = 0.001 a.u., the He~^ (2s) and He~^ (2p) states are 
degenerate and coincide with the energy value of A = 2 ionization threshold of free He~^ ion. As 
A increases, the He^ (2s) and He~^ (2p) states become non-degenerate. Initially, the 2s level 
of He~^ lies energetically below the 2p level for A up to 0.5 a.u. At A = 1.0 a.u., the 2s state 
moves above the 2p level. These results exhibits that an ^incidental degeneracy' takes place for 
2s and 2p states of He~^ at some value of A between 0.5 and 1.0 a.u. and then a 'level crossing' 
occurs between two states having different symmetry properties. Finally these states become 
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degenerate again for A > 100.0 a.u. The incidental degeneracy for He^{2s) and He~^{2p) states 
occur for A in the range 0.5 < A < 1.0. Such incidental degeneracy and subsequent level crossing 
phenomenon are being noted earlier by Sen et. al. m in case of cage confined hydrogen atom 
and by Bhattacharyya et al. [28] in case of helium-like ions within strongly coupled plasma 
environment. 

2. It is seen from both figures-6 and 7 that all the resonance energies (Ej.) are almost unaltered up 
to A = 0.5 a.u., then decrease rapidly up to A = 20.0 a.u. and ultimately saturates. For low 
values of A (say 0.001 a.u.) the resonance energies are identical with those of the free He atom 
whereas for A = 1000.0 a.u. the resonance energies are equal to those of free He atom plus 0.4 

а. u. {i.e. 2.0X Vo). Thus, for a given depth (Vq) of the finite oscillator potential, the variations 
of energies of the bound states and the resonance states of helium w.r.t. the width of the cavity 
(A) are nearly identical. 

The variation of widths (T) of 2sns and 2pnp {^S^) [n = 2 — 5] resonance states w.r.t. A are given in 
figures-8 and 9 respectively. A closer look at the figures 8 and 9 leads us to the following observations. 

1. In general, it can be argued that the variation of widths shows an oscillatory behavior which are 
more pronounced for the higher excited states. It is worthwhile to mention here that recently 
Chakraborty and Ho [3j also reported such oscillation of resonance width (T) for 2s^ (^S) state 
of QD confined helium atom. This feature clearly indicates a possibility of controlling the 
autoionization lifetime of doubly excited states of two-electron ions by tuning the parameters of 
the confining FO potential representing the quantum dot. 

2 . The variations of widths of 2s^ and 2p^ (^5'®) states with respect to A are exactly opposite in 
nature. For state originating from 2s^ configuration, the autoionization width first decreases 
in the range 0.1 < A < 1.0 and after reaching the minima, it shows a large bump around A ~ 

б. 0 a.u. After that it starts to decrease and finally the autoionization width saturates where it 
becomes equal to that of a free He atom. In contrast, for 2p‘^ state, the autoionization width 
hrst increases for 0.1 < A < 1 and then shows a large dip approximately at the same value of A 
for which the 2s^ state shows the bump. 

3. The values of A corresponding to the largest bump in the values of autoionization widths (F) 
of 2sns states and the lowest dip for 2pnp states shift towards the higher values of the cavity 
width (A) for higher excited states. 


Inside the QD i.e. due to the presence of the surrounding FO potential, the charge distribution of 
the impurity ion gets reoriented which produces the behavioral changes as compared to a free ion. 
The nodes or antinodes of the resonance wavefunction lie at the boundary of the QD cavity and the 
interference caused inside the cavity gives rise to the oscillatory behavior of the resonance widths mm- 
The number of nodes or antinodes of the wavefunction increases for high-lying resonance states and 
the oscillation becomes more prominent. 

4 Conclusion 

Structural properties of He atom conhned in a QD, efficiently modeled by a two-parameter weakly 
confining FO type potential, has been investigated in the framework of Stabilization method using 
explicitly correlated Hylleraas-type basis sets. It has been observed that the structure of the impurity 
ion is a sensitive function of the dot size. For very small values of the cavity width, the system be¬ 
haves almost like a free ion whereas, for very high cavity widths, a constant shift equal to the depth 
of the potential are observed in the bound as well as resonance energies. When the dot size becomes 
comparable to the dimensions of the impurity atom, the effects are more pronounced and many re¬ 
markable behaviors such as increase in ionization potential, oscillations in the widths of two-electron 
resonance states, incidental degeneracy and subsequent level-crossing phenomena for one-electron ions 
are observed. The present work is expected to lead to future investigations on the autoionizing states 
of different angular momenta for QD confined two-electron systems. 
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Figure 1: Stabilization diagram for ^5® states of helium atom under quantum cavity. Width of the 
cavity is set at 4.0 a.u. 



Figure 2: The variation of bound state energy eigenvalues w.r.t. the width of the cavity 
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Figure 4: Enlarged view of the stabilization diagram for states of helium atom under quantum 
cavity in the energy range between -8.0 a.u. to -0.64 a.u. Width of the cavity is set at 4.0 a.u. 






Energy (a.u.) 

Figure 5: Density of states and fitted lorentzian for cavity width 4.0 a.u. 



Figure 6: The variation of resonance energies (Er) of 2sns [n = 2 — 5] states and corresponding 

2s and 2p threshold energies with the cavity width (A). 


14 























Figure 7: The variation of resonance energies (Er) of 2pnp [n = 2 — 5] states and corresponding 

2s and 2p threshold energies with the cavity width (A). 
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Figure 8: The variation of resonance width (F) of 2sns [n 
(A). 
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Figure 9: The variation of resonance width (F) of 2pnp [n 
(A). 


2 — 5] states with the cavity width 
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